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Abstract 

One of the strongest constraints on the existence of large, compact, "gravity- 
only" dimensions comes from SN1987a. If the rate of energy loss into these 
putative extra dimensions is too high, then the neutrino pulse from the su- 
pernova will differ from that actually seen. The dominant mechanism for 
the production of Kaluza-Klein gravitons and dilatons in the supernova is 
via gravistrahlung and dilastrahlung from the nucleon-nucleon system. In 
this paper we compute the rates for these processes in a model-independent 
way using low-energy theorems which relate the emissivities to the measured 
nucleon-nucleon cross section. This is possible because for soft gravitons and 
dilatons the leading contribution to the energy-loss rate is from graphs in 
which the gravitational radiation is produced from external nucleon legs. Pre- 
vious calculations neglected these mechanisms. We re-evaluate the bounds on 
toroidally-compactified "gravity-only" dimensions (GODs), and find that con- 
sistency with the observed SN1987a neutrino signal requires that if there are 
two such dimensions then their radius must be less than 1 micron. 
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I. INTRODUCTION 



The conjecture that the four- dimensional universe we perceive is merely a projection 
of a higher-dimensional space dates back at least to the work of Kaluza and Klein in the 
1930s. Such ideas imply that we live in a world in which there are four infinite dimensions, 
and n compactified "extra" dimensions. If this is the case there is the possibility of observ- 
ing Kaluza-Klein (KK) modes, which correspond to excitations of ordinary standard-model 
particles in these extra dimensions. The absence of any such effects in collider experiments 
performed to date then places stringent limits on the size of such extra dimensions, forcing 
them to have inverse radii larger than 1 TeV [0. Until ten years ago the conventional wis- 
dom was that the size of these extra dimensions would be set by the Planck mass, and so 
would be much smaller than this. 

In the last decade this conventional wisdom has been under attack from two different, 
but connected, directions. Firstly, it has been suggested that string theory could lead to 
compactified dimensions which have radii much larger than 1/Mp Secondly, an alter- 

native picture of the extra dimensions has emerged, in which the standard- mo del fields are 
confined to a four-dimensional "brane" , while gravity propagates in the [n + 4)-dimensional 
"bulk" . This has generated considerable excitement because it seems to provide a natural 
solution to the hierarchy problem. Effectively, gravity is much weaker than the three other 
forces because it is diluted via its propagation in the extra dimensions. One simple scenario 
in which this occurs was suggested by Arkani-Hamed et al. Suppose that there are n 

compact "gravity-only" dimensions (GODs) whose radius R is large compared to the scale 
of fundamental physics in the (n -|- 4)-dimensional world, M(4_|_„). For distances r 3> -R, the 
gravitational force will behave like |0: 



where M(4+„) is the fundamental scale of gravity in the (n + 4)-dimensional world. Matching 
to the Newtonian law of gravity, which is the effective theory for long distances and weak 
fields, we see that 



and so if i? ^ l/M(4+n) then the effective four-dimensional Planck scale Mp ~ l/^fG^ will 
be much larger than the underlying {n + 4)-dimensional scale M(4+„) . 

Of course, at distances r R the inverse-square law begins to break down. Choosing 
M(4+„) ~ 1 TeV implies, for a given n, a definite value of the scale R at which this breakdown 
occurs. In the case n = 1 gravity would be modified on solar-system distance scales, but 
for n = 2 the distance R is already smaller than 1 mm. Thus, while the case of one GOD 
is ruled out experimentally, no direct measurement of gravity has, as yet, been made that 
precludes the existence of more than one GOD 



We do not consider Randall and Sundrum GODs Q. 
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An alternative to measuring the gravitational force at distances of order one millimetre 
is to perform collider physics searches at the scale M(4+„) and look for "physics beyond the 
standard model" [pHTl|. It is not immediately apparent that such searches will reveal any 
phenomenologically-significant effects. The only modification to physics at the scale M(4_|_„) 
is in the physics of gravity, and so the violations of the standard model appearing at that 
scale are, naively, only of gravitational strength. However, at such energies a large number 
of graviton KK modes are excited, even though all standard-model fields continue to behave 
as four- dimensional objects. This compensates for the weak coupling between these KK 
modes and matter fields. Production of KK-gravitons, h, and KK-dilatons, 0, as well as 
their virtual contributions lead to deviations from the predictions of the standard model, as 
demonstrated in Refs. [p|-|TT[]. It is possible that such deviations could be detected at future 



colliders lll2l|T3 



Significant constraints have already been placed on the size of compact GODs from 
precision experiments at existing colliders WM. However, these collider bounds are not the 



most stringent ones on the size of these dimensions. Arkani-Hamed et al. pointed out that 
one of the strongest bounds on this physics comes from SN1987a 0. The argument is very 
similar to that used to bound the axion-nucleon coupling strength |[15|- [1^] . The "standard 
model" of supernovae does an exceptionally good job of predicting the duration and shape 
of the neutrino pulse from SN1987a. Any mechanism which leads to significant energy- 
loss from the core of the supernova immediately after bounce will produce a very different 
neutrino-pulse shape, and so will destroy this agreement — as demonstrated explicitly in 



the axion case by Burrows, Brinkmann, and Turner [|1^. Raffelt has proposed a simple 
analytic criterion based on detailed supernova simulations : if any energy- loss mechanism 
has an emissivity greater than 10^^ ergs/g/s then it will remove sufficient energy from the 
explosion to invalidate the current understanding of Type-II supernovae's neutrino signal. 
In the case of interest here the bremsstrahlung processes A^A^ — * NNh (KK-gravistrahlung) 
and A^A^ NNcj) (KK-dilastrahlung), provide an energy-loss mechanism for the supernova. 
Energy generated during the collapse is radiated away into the extra-dimensional degrees 
of freedom, never to return. Demanding that such reactions do not result in the loss of too 
much energy from the supernova allows bounds to be placed on the size of the toroidally- 
compactified GODs. 



Two attempts to place such a constraint have already appeared in the literature |T9|,pO 



Both assumed that the A^A^ scattering amplitude was dominated by a single virtual pion 
exchange — an uncontrolled approximation which we refer to hereafter as "the one-pion- 
exchange approximation" . Although this approximation appears reasonable, since one-pion 
exchange is the long-distance component of the nuclear force, in fact the nuclear physics of 
this problem is much more complicated. Most obviously, the very existence of nuclei tells 



us that the strong nuclear force is non-perturbative. Yet, the calculations of Refs. ||19|J20 
are only valid if one-pion exchange can be treated in Born approximation. Hence, the A^A^ 
interactions which occur during these radiative processes should, instead, be described by 
the full A^A^ t-matrix. 

One might think that such a calculation is model-dependent and complicated. However, 
in the limit that the emitted radiation is soft compared to the incoming A^A^ energy, the 
calculation simplifies dramatically. In this limit there are low-energy theorems for the emis- 
sion of radiation which are analogous to those for photon bremsstrahlung [pi|,p^]. These 
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theorems relate the amplitude for the bremsstrahlung process to the on-shell A^A^ ampli- 
tude. It thus becomes possible to use A^A^ scattering data to calculate low-energy emission 
processes directly. Recently, three of the present authors used theorems of this sort for axial 
bremsstrahlung to perform a more reliable calculation of axion and neutrino bremsstrahlung 
in a nucleon gas ||23|. The resultant rates are model- independent, and rely only on two as- 
sumptions: that the radiation is soft, and that the emission rate from the nuclear gas is 
dominated by two-body collisions. The use of such an approach therefore essentially elimi- 
nates the two-body dynamics as a source of uncertainty in the problem. 

In this paper we will derive model-independent low-energy theorems which relate the 
amplitudes for A^A^ — > NNh and A^A^ — > NN(f) to the on-shell A^A^ scattering amplitude. It 
turns out that at both leading and next-to-leading order (LO and NLO) in the soft expansion 
the cross sections for these production processes can be expressed directly in terms of the 
A^A^ differential cross section. We use the NLO production cross sections to compute the 
energy loss from a supernova due to KK-graviton and KK-dilaton radiation from the two- 
nucleon system. The "Raffelt criterion" explained above then leads directly to a constraint 
on the size of the extra dimensions. 

The paper proceeds as follows. In Section |I| we derive the amplitude for the emission of 
soft gravitational radiation from the A^A^ system. We include all terms which contribute up 
to the first two orders in the soft expansion. We then take the amplitude at NLO in the soft 
expansion and use it to compute the energy loss due to the three processes A^A^ — > NNh, 
NN NN(j), and A^A^ NNg, where g is the four-dimensional graviton. In Section |T] 
we derive formulae, meant for numerical evaluation, which yield the emissivity of a gas of 
nucleons due to these production mechanisms, give our "best" result for the bound on the 
size of the extra dimension^, and discuss the validity of the soft approximation. Section |V 
then examines the limiting cases of non-degenerate and degenerate neutron matter. In these 
limits analytic formulae for the emissivity in terms of a single A^A^ cross section may be 
obtained. We also exhibit the analytic formulae found if one assumes that effective-range 
theory describes the A^A^ cross section, and present the radius bounds given by these analytic 
calculations. Finally, Section ^ presents the results for the supernova emissivity obtained 
with all of these different approaches, and compares and contrasts them with the work of 
Refs. IIJII. 

Many-body effects, although undoubtedly important, were not the focus of Ref. [^, and 
will not be part of our discussion here. While the radiation from the two-body system can 
be computed in a model-independent way, the evaluation of many-body effects is necessarily 
model-dependent, since it relies on selective re-summations of classes of diagrams in the 
many-body problem. One set of many-body effects which can modify the vacuum rates 



for KK-graviton emission was considered in Ref. [0]. We look forward to future work on 
how many-body physics affects the reactions A^A^ NNh and A^A^ — > NN(f), but we have 
not focused on these issues here, since our goal was to provide a firm, model-independent, 
foundation on which future calculations of many-body corrections to these processes can 
build. 



^In general what is actually constrained is the nth root of the product R1R2 ■ ■ ■ Rn, where Ri is 
the radius of the ith extra dimension. 
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II. NN NNh AND NN NN(f) 



A. Graviton and dilaton dynamics in n + 4 dimensions 



For simplicity, we let the n extra dimensions form an n-torus of radius R = L/(27r) 
in each direction. The zero-modes of the (n + 4)-dimensional graviton correspond to the 
massless graviton, h^^^, the {n — 1) vector bosons (which decouple from the matter fields) 
and n{n — l)/2 scalar fields, (p^^. The graviton and scalar field mode expansions have the 
form 



h'^'^ixj) = ^ /i'^'^'^'(x) exp 



.27rj ■ y" 
' L 



[X) exp t- 



2 7Tj ■ y ' 
L 



(3) 



where x are the 4-dimensional coordinates, y are the n-dimensional coordinates, and j is 
an n-dimensional vector representing the momentum of the mode in the extra dimensions. 
Although these modes are massless in the (n -|- 4)-dimensional world this extra-dimensional 
momentum results in an effective mass on the four-dimensional brane, mj, withQ 



L2 



(4) 



We use KK-graviton fields with polarization tensors ejj^i^{k,X), where the extra- 
dimensional momentum is j, while the graviton is propagating in the direction k in our 
world, and has polarization A. Then 



A=l 



1 



(5) 



with ||,|T0| 



2 

mi 



~l~ o VfJ-f ~l nkuki, 

6 \ mi 



(6) 



The Lagrange density that describes the four-dimensional interactions between the KK- 
graviton and KK-dilaton, and matter fields is then P,p!9| 



E 



3{n + 2y ^ 



(7) 



^Note that this predicts the existence of a massless dilaton. This is unphysical, since any mech- 
anism that stabilizes the GODs at size R will provide a mass to the field However, as long 
as this mass is much smaller than the temperature of the nucleon gas in the supernova its precise 
value will not affect our calculation. 
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where T^jy is the conserved energy-momentum tensor of the matter on the 4- dimensional 
brane, and the dilaton field is = 0aa- We have also defined k = ^/32'^■GN, with G^v being 
Newton's gravitational constant^. 

B. Radiation from the NN system in the soft limit 

To compute gravitational radiation from the A^A^ system we must calculate the energy- 
momentum tensor of free nucleons. Neglecting isospin violation in the strong interaction, 
and electroweak effects 

t;^ = - ^ J 7.^ - 9,uC ^ {'n, - 1,) - ^ + ... , (8) 

where the ellipses denote spin-dependent interactions that do not contribute if the momen- 
tum of the emitted radiation is small, and terms that can be made to vanish by the nucleon 
equations of motion. M is the isospin-averaged nucleon mass. This, together with Eq. ([7|), 
defines the leading interaction between a single KK-graviton or KK-dilaton and a nucleon. 

Now we fix our attention on the case of KK-graviton emission, and examine the limit 
when the emitted graviton is soft — which is to say its energy, uj, is much less than all other 
scales in the problem. This means that we certainly want u -C p'^/M, the incoming energy of 
the A^A^ system, although other restrictions will also emerge, as we shall see below. In what 
follows we shall treat uM/p^ as a small, dimensionless, expansion parameter, and compute 
the contribution of the first two orders in this expansion to the amplitude for A^A^ NNh. 

In the limit uM/p^ -C 1 the amplitude for gravistrahlung is dominated by graphs (a)- 
(d) in Fig. In Fig. |T] the incoming nucleon momenta are pi, the outgoing nucleon 
momenta are p'^, p^-, and the KK-graviton four- momentum is denoted by k. We must also 
define Mandelstam variables Sj and sj for the A^A^ interaction, depending on whether the 
nucleons scatter before or after the graviton emission. Likewise, two different Mandelstam 
fs enter, depending on whether the KK-graviton is emitted from nucleon 1 or nucleon 2. 
Denoting 

S^ = {Pl+P2f] Sf = {p[+p'^f (9) 

ti = {pi-p[?; h = {p',-p2?. (10) 

The amplitudes for KK-graviton emission from the four graphs (a)-(d) may then be written 
as: 

^^(^^ = (^^^ (2p[ + k)^{2p[ + k)MPi)u{P2) '^_^ {-t)A{s,M)u{p-'i)u{p-'2), (11) 
^M^} = (^) (V2 + k),i2p', + k)MPi)m') '^_^ i-t)Ais,,t,)u{p-l)u{p-'2), (12) 
iM^^J = (^) HprMp-,'){-^)A{sJ,t,)^—j^{2p, - k),{2p^ - k) Mp{)u{p-',\ (13) 

^Our K is a factor \/2 larger than that employed in Ref. |^] , and our graviton polarization tensor 
is a factor of 2 smaller. This change does not affect the answer computed for any observable. 
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(«) A A 
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Pi ^ P{ Pi P{ 

(b) A A (^) 

N S 
S \ 

h,(f> h,(f> 



— Pi 

k 

h,(t> 

— Pi 

FIG. 1. The leading diagrams contributing to processes NN — > NNh and NN NN(j). Nucle- 
ons are denoted by solid lines and the KK-modes h or cp are denoted by dashed lines. Solid squares 
denote an insertion of the single-nucleon energy-momentum tensor, while solid ovals containing A 
denote an insertion of the full NN scattering amplitude. The solid oval containing T denotes the 
non-pole vertex required for the sum of diagrams to satisfy d^M^'^ = 0. 




^MjfJ = (^) u{p-[')u{p2'){-t)A{sf, h) ^ _ ^ (2P2 - k),{2p2 - kUiPiMP2)- (14) 

This amplitude begins at 0(x^^), with the small parameter x is meant as a mnemonic for 
any of the quantities \k\, uj, and mj. The u's are Dirac spinors, normalized so that uu = 1. 

Note that here the amplitude A is the on-shell A^A^ amplitude projected onto positive- 
energy A^A^ states. One might think that is necessary to include in the expressions (|TT|)-(|T^ 
the fact that one leg in the A^A^ interaction is off its mass shell. However, it is clear that the 
difference between using the on-shell A^A^ amplitude and the one that is "correct" for the 
kinematics of Fig. is an 0{x^) effect. The low-energy theorem of Low |21|] for A^A^ — >■ NN'j, 



and that of Adler and Dothan p2[ for axial bremsstrahlung, are both based upon the idea 
that we can construct the leading-order bremsstrahlung amplitude using the on-shell A^A^ 
amplitude, and then determine the terms of 0{x^) using current conservation. As we shall 
see below, a similar low-energy theorem exists in this case, and so we need not bother to 
include the off-mass-shell dependence of the A^A^ amplitude in Eqs. ([TT|)-([1^. In fact, in 
general the off-shell dependence of the A^A^ amplitude is the result of a particular choice 
for the nucleon interpolating field, and so it can always be changed by field redefinitions. 
S-matrix elements are unaffected by such field transformations 

We now simplify the above expression, by using the fact that A is projected onto positive- 
energy states, and retaining only the spin-independent terms up to 0(x°)- Although it is 
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possible to keep track of the spin- dependent terms, we will drop them here, since they always 
appear at 0(x°), and so, after summing over nucleon spins, do not make a contribution to 
the production cross section until 0(x°)- This is two powers of x beyond the leading term 
in the cross section. Thus, excluding these terms, the result is: 



2p\ - k + k"^ 



2M2 



k-pi 



2p'2 -k + k^ 
1 



2M2 



i2p,-k)^i2p,-k). 



2pi ■ k — k"^ 

where we have also used the fact that the external nucleon legs are on-shell. 



(15) 
(16) 
(17) 
(18) 



C. Conserved M, 



Observe that if we sum graphs (a)-(d) to get M^^, then the result is not a conserved quan- 



tity. To remedy this, we generalize the prescription used in photon bremsstrahlung pi| , |22 
writing: 



k^'M, 



-k^T^, + 0{x^ 



and then construct 



(19) 



(20) 



By the same arguments that apply in the case of the Low ||2T| and Adler-Dothan 



theorems, M^J will be the correct amphtude up to spin-dependent pieces of order x° and 
spin-independent pieces of O(x^) Q- 

A brief calculation then reveals that the spin-independent piece of the amplitude for 
NN NNh is, up to 0(x°): 



= - 



A{Sf,t2 



1 



2pi ■ k — k"^ 



{2p,-k)i2p,-k)^ + A{sfM 



1 



2p2 ■ k — k"^ 



{2p,-k)i2p, 



{2p\ + k)^ {2p[ + k)^ ^ , ] , ^ Ms„ t^) - {2p', + k)^ {2p'^ + k),— ^ 



2p'i -k + k"^ 
+ 2r]^,^A{s, t) + 4(pi + P2)m(Pi + P2) 



2p'2 -k + k^ 



A{si,ti 



dA{s, t) 



ds 



+ 4(P2 -P2)/.(P2 -P2)i. 



dAis, t) 



dt 



t=t. 



^Structures involving the spin-matrix a^j^u which are not constrained by current conservation can 
contribute to M^i, at order x^ . However, the effect of such terms on the observables computed 
here is suppressed by two powers of x compared to the leading soft result. 
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The classical amplitude for the emission of gravitational radiation from a system of 
colliding particles was derived many years ago by Weinberg |26|. We may now recover this 
result by retaining only the 0{x~^) parts of Eq. (|21|): 

Mjr.'^ = f E vm,p^.-^As, t), (22) 

where j, the particle index, runs from 1 to 4, with j = 1, 2 for incoming nucleons and j = 3, 4 
for outgoing nucleons, and 

"^'-{-1 ifj = 3,4; ^''"^^ 
t^^. (24) 

Weinberg's derivation of Eq. (|22|) is entirely classical, which is not surprising, since the 
graphs that contribute to the amplitude only involve radiation from external legs, and so they 
correspond to classical bremsstrahlung from the gravitational charges accelerated during the 
NN collision. It is also worth pointing out that the Weinberg result is only valid if the A^A^ 
t-matrix is slowly varying over the energy region of interest^. 

Next observe that the conservation of M — k^M^^^j = — implies the following relations: 

y 

Mo. = Mji, (25) 

Moo = -^Mji. (26) 

This obviates the necessity to actually calculate these "0" components. These relations may 
be formalized by defining the symbol Wfj,uij as follows: 



W^,yij = rii^rijy - rii^riou— - Vot^Vjiy— + Vo^iVou^^- (27) 

U! UJ UJ 



If M^j, is conserved then we need only calculate the space-space components of this tensor. 
Other portions can be generated via the relation: 

M^, = W^,,,M'^. (28) 

In particular, if we work in the centre-of-mass frame, then at leading order in the nucleon 
velocity the space-space components of M^^ have a very simple form: 

Mt'^=Ais, i) [p,p, - p[p'^] . (29) 



^Weinberg's basic result actually omits the amplitude A. He assumes the scattering that leads to 
the emission of the radiation is isotropic, and then includes anisotropy in the two-body scattering 
by multiplying by the two-particle differential cross section. 
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This result differs from Eq. (^) by relativistic corrections, which are suppressed by two 
powers of p/M. In addition, the next-to-leading order (NLO) amplitude (21) contains cor- 
rections to Eq. ( [29| ) which are suppressed by u /M. Since uj <^ p^/M, by assumption, the 
expansion in p/M is less convergent than the one in uj/M . 

If we continue to neglect these effects of order p^/M^ and uj/M, then we can also write 
down a compact formula for the next-order correction to Mij-. 



(0) 



--5ijA{s, t) + {pi- p[){pj - p.] 



dA 



s,t) 



t=t. 



(30) 



where, once again, we have chosen to work in the centre-of-mass frame. 

Note that working only at leading order involves neglecting the pieces of the amplitude 
M^j, corresponding to variations of the NN t-matrix in the Mandelstam variables s and 
t. This is clearly not a good approximation at low energies, where the A^A^ t-matrix varies 
rapidly. At higher energies the scale of variation of A can be thought of as m^, the pion mass, 
and so we would also like to have uo -C m.,^. In fact, this is probably a rather conservative 
estimate, since empirically the energy and angle- dependence of the A^A^ cross section is 
quite weak at these higher energies. However, a proper test of the convergence of the soft 
expansion is really required. 

One obvious way to perform such a test is to compute the LO result (2£) and then 
compare it to the result at NLO. In Section [II E| we will show that in fact the NLO correction 
vanishes, and so the LO amplitude (E^) is more reliable than one might naively expect. 



D. Result for energy loss from the NN system to leading order 

In this Section we will compute the energy-loss for a given graviton KK-mode j. In 
Section |ITI| we sum over all modes in order to compute the total rate of energy emission from 
a gas of nucleons. 

— * 

The energy loss due to graviton emission into a given KK mode j is given by: 

dei = ^J^^ E [M'-'ei.ik, A)J ej,^{k, X)M^\ (31) 

(Strictly speaking this is not an energy loss, since we need to multiply by the density of 
initial and final A^A^ states, however, throughout this section we employ this somewhat lax 
terminology.) 

Using the symbol W, defined in Eq. (|27|) , we may rewrite this as 



di = ^^^M^^W^-^bI^.W-^'^^M,,, (32) 

where B is defined by Eq. (^). Since k^W^^ki = by construction, only the pieces of B 
proportional to 7] symbols survive. After some calculation we see that: 



-;^^,uijB^'"^'^Wal3kl = - [VikVjl + VilVjk - ^VijVkl 

(33) 
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At this point we may integrate over the final-state graviton direction k, since at LO Mij 
does not depend on k. This gives: 



1 f dk 19 1 



rrA /\\ 14 \ 2 

+ {-^{^ik^ji + ^ii^jk) - -^^ijhij + "^^(^^'^^J' + ^^'^i*^ + ^ijhi)- (34) 

Using the LO M^^, Eq. (p9D, we finally obtain the following energy loss per unit frequency 
interval, for radiation into a KK mode of given mass mj, 

den k8Gr,4f\ , .,2 ^9 11 "^j 2"^| 
aci; cij ovr iW^ \ 18 9 cj^ 9 cj^ ' 

where 6cm is the centre-of-mass scattering angle: cos 6*0111 = P ■ P', and = {p^ + p'^)/2. 

This argument can also be pursued for 4-dimensional gravitons, provided one modifies the 
polarization tensor B appropriately. The result for energy-loss due to 4D-graviton emission 
is: 

2 

dsg %Gn^P^ .2a 1/112 (oa\ 



as derived by Weinberg ||2^. Note that this is not the mj — > limit of Eq. (pSf), since the 
4D-graviton has fewer degrees of freedom, due to its 4D-masslessness. Note also that here, 
and in Eqs. (^Sj) and (|38|), we have made the replacement 

P^+/-2pV'cOs2^em 2^%in2^em, (37) 



in order to make contact with Ref. [|26| and simplify the formulae. This corresponds to 
evaluating the energy of the NN system at the point s. The error introduced by this 
replacement is of order M^cj^, and so is suppressed by two powers of x- 

Meanwhile, the result for energy loss due to KK-dilaton emission is also obtained by an 
analogous argument. It is: 



dsA, k^G^Ap^ 



■2 ^ / ^2 



2 



sin^^eml^pTTT— ^ • (38) 



duj uj Stt M2 ' 9(n + 2) 



E. Next-to-leading order calculation 

These energy- loss rates represent the leading-order result, which we shall denote by 
0(x~^)- To calculate the correction to this energy- loss rate which arises at NLO we must 
contract the amplitude (^) with the correction (0). Since — at least in the non-relativistic 
limit — the space-space components of the NLO amplitude do not depend on the direction 
of the emitted radiation, we may use the steps (pl|)-(|3^) employed above to show that the 
piece of ^ which arises at 0(x~^) vanishes. 
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To do this we first note that the trace of M*^ over its spatial components, Mj;^ ^\ is 
0{x^). Therefore, since we are trying to compute pieces of ^ that are 0{x~^) we can 

neglect terms such as MaMjj. It also follows that the contraction of the 6ij piece of m/°^ 

with M^;"^"* only produces pieces of the energy-loss rate that are of higher order than we are 
concerned with here. 

A brief calculation then shows that the contraction of the piece of M^j'^ proportional to 

^ with also only yields pieces of 0(x°). Thus, it follows that the first corrections to 

the expressions (|35D, (p6D, and ( |38D are suppressed by two powers of the small parameter x, 
which we can now identify more precisely as 

X = ^- (39) 

It is worthwhile stressing that the NLO correction to these results would not be zero if we 
had used, say, the initial-state energy p'^/M in these formulae instead of the average energy 
of the A^A^ state: p'^/M. It must also be remembered that at 0(x°) there exist contributions 
to the energy-loss rates which are not constrained by the dual requirements of the on-shell 
A^A^ data and the conservation of M^^,. 



III. GRAVITATIONAL EMISSIVITY OF A GAS OF NUCLEONS 

A. Non-Relativistic Nucleons at Arbitrary Degeneracy 

The gravitational emissivity of a nucleon gas due to the two-body reaction a + 6 — *■ 
a + b + X , where the labels a and b can represent either a neutron or a proton, and X 
is the type of gravitational radiation (massless 4-D gravitons (g) or massive Kaluza-Klein 
gravitons (h) or dilatons (0)), is given by the general formula 



dSx 
dt 



■ -j-j. d^pi d^p[ 



S UE,)h{E,){l - fa{E[)){l - ME',)) 



dp 

{27r)'6\p,+p,-p[~p',-k)^, (40) 

Here is the energy-loss per unit frequency due to A^A^ — > NNX computed above, 
summed over initial and final spins, for the emission of KK-gravitons [see Eq. (|35D], ordinary 
gravitons [see Eq. (^6l)], or KK-dilatons [see Eq. (^)]. Again, as noted above, Eqs. (|35D - 
(^) are not really energy-loss rates, since they do not include the A^A^ density of states. 
However, when integrated as in Eq. (^0|) they will yield the emissivity. Meanwhile, fa and fb 
are the nucleon Fermi distribution functions for the species a and b. The symmetry factor 
S = [1/(1 + Sab)]'^ is different in the np and nn cases. The sum over KK modes clearly does 
not apply to the case when X is the four- dimensional graviton. 

The formula (|40|) applies equally well to the np and nn cases. Although we could proceed 
in complete generality and consider a nucleon gas with different densities of neutrons and 
protons, we will specialize our treatment at this point to the nn case, since the formulae 
are more transparent there. The proton fraction is generically small in supernovae, and so 
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the inclusion of np scattering — which is easily done using steps analogous to those discussed 
below — should make little difference to the radius bound we obtain. 

In the soft-radiation limit we may neglect the momentum k in the three-momentum 
conserving delta-function of Eq. (|40|) . Spherical symmetry and energy-momentum conser- 
vation can also be exploited, thereby eliminating nine of the fifteen integrals appearing in 
Eq. (0). Defining total, initial relative, and final relative momenta P, p, and p' we move to 
dimensionless variables [ITB 



P' P' , P" .... 
Up = ; Ur = ; = , 4i 

8MT' 2MT' 2MT' ^ ^ 

where T is the temperature of the neutron gas. 

Denoting cos(6') = p- P, cos{9') = p' ■ P and (p to be the angle between p and p' projected 
onto the P plane, we are left with a six-dimensional integral: 

ur ~ ^ 57r6 J 5^ 



S ^— ^ I dur I d(cos9) J dup j ^ du[. j d(cos9') 



(-277 



Uy'u](\'/'U' exlVK -<)] - - /2) /„ ^ Sin^^em \Ai9,m, 2Tu)\' . 







(42) 

Here fi = [exp(Mj — y) + with ui^2 = up + Ur ± 2^upUr cos 9, and y = fi/T. The f-s 
are defined similarly in the final state, with u[ 2 = mJ. + ± 2y/ u[.up cos 9' . We have also 
defined 5j = mj/ (2T) and u = {ur + u'j.)/2, and written ^ as a function of the centre-of-mass 
scattering angle and the centre-of-mass energy. An implicit summation of over both 
initial and final nucleon spins is to be understood here, and throughout what follows. The 
function ^ is different depending on the species X under consideration: 

e.N = 1, (43) 
aw = v/r^(^ + ^x^ + ^x^), (44) 

2\5/2 



f*w = 9(;rT2)('-^) • 

The dependence of the integrand on the angle (p enters solely through the centre-of-mass 
scattering angle ^cm, which is related to 9, 9' and (p via 

cos 9cm = cos 9 cos 9' + sin 9 sin 9' cos . (46) 

Now we wish to perform the sum over KK-graviton modes. If the modes are narrowly 
spaced we may rewrite: 

— > P^^n j dmmJ'-^, (47) 



where fi„ is the surface area of the n-dimensional sphere. The order of integration may then 
be interchanged, and the integral of the function / over the graviton masses evaluated. We 
may also change variables to u'^ and a = 00 /T, so that ultimately we obtain: 
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d£x _ 4:Gn 
If ~ 3 



/* OO P OO 

{2TY^/\TRYTi^M^'^hn da du',K^''\a,u',), (48) 

Jq Jo 



where the dimensionless kernel K^^'' is 



/ ^sm'9,,^\A{9,^,T[2u[. + a/2])\' / dup ^f^l - f[){l - f2), (49) 
Jo In Jo 

and 

_3n^ + 24^2 + 55n + 42 1 
" - (n + 2)(ri + 3)(n + 5) T (^) ' ^^^^ 

if we consider the radiation of KK-gravitons and KK-dilatons. Note that the gravitons alone 
give 

_ 3n^ + 18n + 19 1 
"~ (n + 5)(n + 3) r(i^)' ^ ^ 

which, for n = 2, differs from the full result by about 1.5%. 

Note that the number of GODs, n, enters here through the power of the radiated energy 
which appears in the kernel, and through the pre-factor. Different models of the extra- 
dimensional physics would change the density of states in the integral over the KK-graviton 
modes, and so would give a different overall pre-factor, but the structure of the answer for the 
emissivity would be essentially the same. The crucial fact about the extra dimensions which 
makes the emissivity (|48| ) sizable during SN1987a is that for i? ~ 1 mm the dimensionless 
quantity TR is of order 10^^, and so many KK modes are excited during the A^A^ collisions 
that take place in the violent supernova environment. 



B. Numerical results and bounds on the size of extra dimensions 

During the supernova, the newly-born neutron star (proto-neutron star) has a typical 
temperature of 40-60 MeV and a high baryon density corresponding to a neutron chemical 
potential fin — 50 — 100 MeV in its inner regions. Thus, matter has f^n/T ~ 1, and so is 
neither degenerate nor non-degenerate. Under such conditions numerical evaluation of the 
phase-space integrals in Eq. (H) is necessary. 

This is done using an NN amplitude reconstructed from the experimental A^A^ phase 
shifts [^. Although there is data on nn scattering only at very low energies, the wealth 



of experimental pp data may be used to infer T = 1 phase shifts. By isospin symmetry, 
these phase shifts should be the same as those for nn scattering, up to small corrections for 
charge-independence breaking. Once a partial-wave decomposition is performed on A the 
(f) integral can be performed analytically, leaving us to evaluate a five- dimensional integral 
numerically. The results for the emissivity in neutron matter at nuclear matter density, 
= ^0 = 0.16 fm~^, in the case of two GODs of radius 1 mm are represented by the solid 
line in Fig. |[ 



14 



10" 



10" 



£10" 

DC 



10" 



10" 



\ 








n=2 




\ 






\\ 






\\ 






\\ 






■ \\ 




- 


V\ 






V\ 






\\ 






7 \ 






\ ~ \ 






\^ \ 




"V. 




10 20 30 40 50 10 20 30 40 50 
T [MeV] T [MeV] 

FIG. 2. Bounds on the radius of the extra dimensions as a function of temperature for neutron 
matter at = 0.16 fm~^ (solid line) and nj^ = 0.48 fm~^ (dot-dashed line). The left panel gives 
the results for n = 2 and the right panel those for n = 3. 



With these numerical results in hand we use the simplistic upper bound on the power per 
unit volume lost from SN1987a of £nn = 10^^ ergs/g/s [l^ to place a bound on the radius 
R of the extra dimensions. The resultant bounds are shown, as a function of temperature 
for the cases n = 2 and n = 3, in Fig. |^. The solid line is the result at = 0.16 fm~^, 
while the dot-dashed line is the result at three times nuclear matter density. 

Assuming that our analysis applies to SN1987a, and using T = 30 MeV and a neutron 
number density equal to the density of nuclear matter, our "best" bound on the compacti- 
fication radius is: 

i? < 7.1 X 10""^ mm for n = 2, 

R<8.5x 10"^ mm for n = 3. (52) 
These bounds are a factor of two to three weaker than those obtained by CuUen and Perel- 



stein, who found |T9 



i? < 3 X 10 ^ mm for n = 2, 

R<Ax 10~^ mm for n = 3. (53) 



The source of the difference between the bounds (q^) and our bound ( |52D will be discussed 
in Section |V|. 

It is interesting to examine the energy distribution of both the radiated particles and that 
of the emitting nucleons, in order to determine if the soft-radiation approximation is valid. 
To do this, note that the kernel K^^-^a, u'^,), defined in Eq. (^9]), may be integrated over u'^, or 
a in order to generate a distribution in the radiation-energy variable or the nucleon-energy 
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energy [MeV] energy [MeV] energy [MeV] 

FIG. 3. Two different partial integrations of the dimensionless kernel, K^^> f ^ \ , for the 
case n = 2, and nuclear matter density. The solid curve corresponds to the final-state nn energy 
distribution, while the dot-dashed curve corresponds to the radiated particle energy distribution. 
For T = 10 MeV and T = 30 MeV the curves have been multiplied by a factor of 10 and 100, 
respectively. 

variable. (Recall that a = u/T and u'^. = E^^f/ (2T).) The results of this procedure at three 
different temperatures are displayed in Fig. as a function of the energy of the radiation 
and the final-state energy of the nucleons. 

Figure ^ shows that at low temperature the typical radiated energy is much less than the 
typical total energy of the final-state neutrons, and so soft-radiation theorems should work 
very well. This is because in the degenerate limit the nucleon energy is set by p^?, while 
the radiation's energy is proportional to T. However, as the temperature is increased, and 
the neutrons become non-degenerate, both energies scale with T. At temperatures relevant 
to SN1987a the radiation is seen to carry off energy comparable to that of the final-state 
nucleons. Note that if the neutron matter were denser then the Fermi momentum of the 
nucleons would be higher, and so the temperature where the peaks of the energy distributions 
coincide would be somewhat larger. On the other hand, as the number of extra dimensions 
rises the power of a appearing in the kernel K*^") increases, and so the radiation would tend 
to carry away even more of the NN state's energy. 

The right-most panel of Fig. |^ represents the conditions used in deriving the radius 
bounds (|52D . We infer that roughly half of the initial energy of a typical A^A^ pair is 
radiated away to the GODs. In other words, if gravitational radiation was emitted into 
extra dimensions during SN1987a it may not have been soft. The parameter x, defined in 
Eq. (^), which we expect will control the expansion is of order 2/3 for the conditions of 
interest. This is a rather large value for our "small" parameter. However, the first non- 
vanishing correction to the emissivity is suppressed by two powers of x, which is somewhat 
reassuring. 
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IV. LIMITING CASES AND ANALYTIC FORMULAE 



While the core region of SN1987a is neither degenerate nor non-degenerate neutron 
matter |]l6l , it is illuminating to consider dilastrahlung and gravistrahlung processes in these 



extreme cases. As we shall show in this section, in both limits analytic formulae for the 
emissivity in terms of one A^A^ cross section can be obtained, provided certain observations 
about the differential cross section are approximately true. Alternatively, effective-range 
theory can be used to provide a form for ^ which can be integrated analytically. In a 
supernova the neutron momentum distribution lies essentially mid-way between degenerate 
and non-degenerate, and so these limiting cases provide bounds on the true emissivity. 

A. A non-degenerate neutron gas 

To analyze the non-degenerate limit of neutron matter the final-state blocking factors 
(1 — /j'), are set equal to unity. In other words, all of the final-state phase space is open 
to the scattered neutrons. In addition, the initial-state Fermi distributions are replaced by 
Maxwell-Boltzmann distributions. These simplifications allow the emissivity to be written 
in terms of integrals over the initial centre-of-mass energy distribution, the energy of the 
KK-dilaton or KK-graviton, and the scattering angle 



dSnn 8Gn T'^ 2 /nrTi\nt. n+i r(?2 + 1 



dt 3 VM'-^ r n + ^ 



2 

/■oo fl firr 

/ dTT^+'e-^ rf(cos^em) sin^^,^— — (Tr,^,^). (54) 

Jo J -I Uiicm 

The variable r = E^^^^^/T is the ratio of the initial nn centre-of-mass kinetic energy to the 
temperature (r = 2ur), and n^ is the number density of neutrons. We have also used the 
relationship between A and the A^A^ differential cross section: 

da M^vrVl' 



n'^ + 9n + 23 
(2n + 7)(2n + 5) 



dn (27r)f 



(55) 



In order to get the form (0) we had to replace the energy at which A is evaluated, 2Tu, 
by 2Tu, the initial-state energy. As we shall see below, this replacement makes very little 
difference to the final result for Snn- In contrast, if we had employed the initial-state A^A^ 
energy p'^/M, instead of the average energy p'^/M, in the pre-factor for the energy- loss rates 
then the emissivity would have been larger by roughly a factor of two. However, as 



discussed in Sections |II D| and [II E| , the choice p^ /M is the one that pushes the error in the 
LO result to 0{x^). Therefore, the formula ( |5^ is an excellent approximation to the NLO 
result for the emissivity in the non-degenerate limit. 

We now consider two analytic forms for the nn scattering cross section. 

1. Nondegenemte neutron gas: energy-independent ^ 

The nn differential cross section is, in fact, essentially energy-independent in the region of 
interest. The angular dependence is also weak, except at small angles, and the contribution 
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of these angular regions is suppressed by the presence of sin^6'cm- The emissivity can thus 
be estimated as 

dSnn _ 16G^ /^.oT^n .2 n„ F + 1) r(n + 5) + 9n + 23 



9 VM^ ^ " ^ ° r(n+|) (2n + 7)(2n + 5)' 



The nn total cross section can be extracted from pp scattering data, provided that Coulomb 
effects are removed. Doing this, and evaluating cTq" at the peak of the function r"'"'"^e~'^, 
leads to a™ = 25 mb, which in turn yields bounds: 

i? < 7.0 X 10~^ mm for n = 2, 

R<8.3x 10"^ mm for n = 3, (57) 



at T = 30 MeV and nN = n^. These limits are very close to the "best" bounds (p2D. 
This indicates that, under the conditions examined here, the neutron gas is very nearly 
non-degenerate, and the neutrons scatter via a cross section which is almost independent of 
angle and energy. Note that even had the nn cross section not been so flat the result ( |57D 
could still be taken as a result of naive dimensional analysis, since any "typical" A^A^ cross 
section will give numbers of this order of magnitude. Its success as such an estimate hinges 
on the fact that Eq. (^) approximates the temperature dependence of the full result very 
well (see Fig. g). 



2. Nondegenerate neutron gas: ^ From Effective-range theory 

Effective-range theory (ER) describes very low-energy nucleon-nucleon scattering by a 
unique momentum expansion about \p\ = |28|. Empirically, the first two terms in the 



expansion of \p\ cot 5 describe the S-wave cross section well, even at energies greater than 
the naive breakdown scale of the expansion, which is \p \ = m-,^/2. The nn non-relativistic 
differential cross section in ER is, neglecting contributions beyond the effective range. 



dVt cm 



i\p 

for 



1^1 > 



1 



(58) 



where ao = —18.5 ± 0.4 fm and ro = 2.80 ± 0.11 fm are the scattering length and effective 
range for nn scattering in the ^5*0 channel As the typical momentum involved in the 
scattering events of interest is much larger than the inverse scattering length we take the 
second form for the differential cross section. Inserting this in Eq. ( p^ leads to a simple 
expression for the emissivity: 



d£. 



ER 

nn 



dt 



1024G7V / n^Nn , itH 2 



r(n + l)r(n + 4) 



9 



N' 



M^rlT (n + 
+ 9n + 23 



(2n + 7)(2n + 5) 



X 



e"r(-(n + 3),x) 



(59) 
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with: 



and r(x, y) the incomplete Gamma function, which obeys r(2;,0) = T {z). The bounds 
obtained from the "Raffeh criterion" £ < 10^^ ergs/g/s at hn = uq and T = 30 MeV are 

R<2Ax 10"^ mm for n = 2, 

i? < 2.1 X 10~*^ mm for n = 3. (61) 

These bounds are much weaker than the bounds (|57|) and (|52|) . As expected, they give 
an upper bound on R, since the cross section (^) includes only S-waves, and other partial 
waves also contribute significantly to Uq" in the energy regime of interest. 



B. A degenerate neutron gas 



Supernovae are, of course, systems far from the degenerate limit. However, the degenerate 



limit is of particular interest for neutron-star cooling, as stressed by Friman and Maxwell ||30 



and Iwamoto |^ in their respective studies of uu and axion radiation. In this section we 



employ some of the techniques used in Refs. fS^JSTI to derive analytic formulae for this case, 
where /i/T ^ 1. 

Radiation from a degenerate neutron gas can only arise from scatterings involving neu- 
trons near the Fermi surface in both the initial and final states. In the soft-radiation limit, 
the momentum of each nucleon involved in the scattering event is constrained to lie on the 
Fermi surface. This approximation leaves only angular and energy integrals in the calcula- 
tion of the emissivity. The energy integrals are those that enter Fermi liquid theory, and 
are well known |^. In particular, the energy integral that occurs in the case of 2 — 2 
scattering is 



/ 



dE, dE2 dE[ dE', hUl - /0(1 - f'MEi + E, + E[ + E',-iu)= ^ ^(^' + 4^'^') , (52) 

e^'-^ — 1 



The u integral can then be done analytically, leaving two angular integrations, and a result, 
for a degenerate neutron gas, of: 



dS„ 



dt 



Gm {RTT K Yn 



r (n + 4) TT^ 
36 

C(n + 2) 



dt ( 



(63) 



(n + 2)(n + 3)' 

where C,[x) is the Riemann Zeta function. The kinematic variables in the centre-of-mass 
frame are related to the integration variables t and 0, defined in the rest frame of the Fermi 

gas by 

2t2 



cos6'r 



\V\-V2\ 



1 - 

2p; 



4p2 



COS( 



(64) 
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where pp is the Fermi momentum, related to the neutron number density hj un = p^/{37r'^). 
Note that this method of calculation is only valid to LO in the soft expansion, since we have 
implicitly assumed that all nucleon momenta are on the Fermi surface. 

Once again, we can now consider two approximate forms for the nn scattering cross 
section. Firstly, evaluation of the integral appearing in Eq. (pSf ) is straightforward if the 
differential cross section is independent of energy and angle. We find that 



dt 



nn 




16 r (n + 4) p\ 

135 ■ 



(65) 



On the other hand, explicit evaluation of Eq. ( |63|) with the differential cross section of 
ER, Eq. (H), yields 



^ = {RTf K Yn 4^'^"- ^ + 4) p| /p^X 



dt 

Dim 



log X + vT+~a? 



X' 



'VT+ 



x^ 



(66) 



We will not use these formulae to place bounds on R from SN1987a. The nucleonic gas 
of SN1987a was far from degenerate, and so such bounds are not particularly meaningful. 



V. DISCUSSION 

In Figure ^ we present numerical results for the KK-graviton emissivity of a neutron gas, 
due to neutron- neutron collisions in the soft-radiation limit, at nuclear matter density. The 
solid curve shows numerical results obtained by using the "measured", angle-dependent, 
elastic neutron-neutron differential cross section and the full Fermi distributions for the 
initial- and final-state nucleons. This numerical result, which is exact at the two-body level, 
provided the radiation is soft, is to be compared with the result given by Eq. (^) when the 
final-state blocking is ignored and the initial-state Fermi-Dirac distributions are replaced by 
Maxwell-Boltzmann distributions. This result, which is valid in the non-degenerate limit, is 
represented by the dotted curve. The dashed line is the analytic result (|56D, that is obtained 
by neglecting the energy and angular dependence of the nn differential cross section. This 
function depends only on the "measured" total nn cross section. Finally, the dot-dashed 



line is the result quoted by CuUen and Perelstein for n = 2 in Ref . |T9[ : 

3.7 X 10^^ ergs/g/s PuT^^^vRmm, (67) 



dt 



where pi4 is the mass density in units of 10^^ g/cm^. Eq. ( |67D was derived using the one- 
pion-exchange approximation for the matrix element and the non-degenerate limit for the 
neutron gas. 

Figure § shows that the conditions examined here are close to the non-degenerate limit. 
Calculations employing this limit produce emissivities only 15% larger than the full calcu- 
lation using Eq. (^8[) . Furthermore, the approximation that the A^A^ cross section is fiat in 
both energy and angle works exceptionally well. The effective-range theory result — which 
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FIG. 4. Emissivity of SN 1987a due to the reactions NN NNcf> and NN NNh, calculated 
assuming there are two GODs of radius R = 1 mm. The solid line is the full calculation, employing 



Eq. (4^), while the dotted and dashed curves represent the different approximations to this discussed 
in the text. The dot- dashed line is the result of Cullen and Perelstein, Ref. 



is not displayed in Fig. ^| — has the wrong temperature dependence. The "full" tempera- 
ture dependence — at least above T ^ 10 MeV — is very close to the result from a fiat cross 
section and the non-degenerate limit: £aT^'^. The emissivity increases with temperature 
more slowly than this if the effective-range theory is employed because the cross section ( |58D 
drops with increasing nn energy. 

The temperature dependence £ a T^'^ is the same as that found in Refs. |ll9|,p0|) although 



Fig. ^ shows that the overall size of the emissivity ( |H7| ) is significantly larger. Refs. |T^.2(] 
both calculated KK-graviton radiation from the NN system by coupling gravitons in all 
possible ways to the one-pion-exchange graph. The resultant diagrams are depicted in 
Fig. |. The authors of both Refs. [|19|,^ had graphs (a)-(d) giving zero, with graphs (e)-(g) 
supplying the full result. This apparently occurs because the square of the matrix element 
was only calculated up to order p^/M^, and graphs (a)-(d) do not contribute until 
If diagrams (e)-(g) are the full amplitude for A^A^ NNh and NN —>■ NNcj) then, in the 
limit m\ <^ MT the matrix element is constant with uj and the NN system's energyQ. 

Our work; here shows that in the soft regime, where uj <^ p^/M, graphs (a)-(d) will 
dominate over graphs (e)-(g). Consequently, if Refs. |TP|,|2D| had examined a degenerate 



^Replacing the full one-pion-exchange-approximation matrix element by its value in the m,r 
limit can be a poor approximation. In the case of axial radiation it enhances the emissivity by as 
much as a factor of two, as discussed in Refs. [^ ,33]. 
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FIG. 5. The seven graphs, which, together with their fermion interchange partners, were cal- 
culated to obtain the KK-graviton and KK-dilaton emissivity in previous works. The solid lines 
are nucleons, the dashed line is the emitted KK-graviton or KK-dilaton, and the dot- dashed line 
represents the pion. 



system, where the gravitational radiation is soft, they would have missed the dominant 
contribution to the emissivity, and so obtained the wrong temperature dependence. On 
the other hand, in the supernova case the matter is far from being degenerate, and as the 
discussion below Fig. |^ shows, the radiation is not especially soft. In fact, naive dimensional 
analysis shows that in the non-degenerate limit: 

So, it should come as no surprise that our graphs, which represent radiation from external 
legs, give the same temperature dependence as graphs (e)-(g) of Fig. ^. Graphs (e)-(g) 
correspond to radiation from internal lines, and appear first at 0{x^) in our approach, but 
they will be accompanied there by many other diagrams. One group of additional graphs at 
0{x^) will represent the gravitational radiation coupling to short- distance components of the 
A^A^ force, while another group involves the inclusion of final and initial-state interactions 
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for the neutrons. Both of these effects should give sizeable corrections to graphs (e)-(g), and 
neither was included in the calculations that led to the result (|67|) . 

Thus, we would argue that the result ( |67D is best viewed as naive dimensional analysis 
of the emissivity £nn due to NN NNh. In this context, we note that either Eq. (|56D 
or the graphs (e)-(g) of Fig. |^ can be used to estimate the emissivity for a non-degenerate 
gas of nucleons, and the dependence they predict on the temperature, and on the size and 
number of GODs, will be correct. The coefficient multiplying the overall result will differ, 
but it will involve "typical nuclear scales" and so it is, perhaps, not surprising that the 
resultant emissivities are within an order of magnitude of each other. Ultimately though. 



our result for Snn is significantly smaller than that of Refs. ||T9|,^, and our bounds on R are 
correspondingly weaker. 



The bounds ( ^2] ) obtained in this work on the type of extra dimensions proposed in 
Refs. 0-0] are model-independent, and are valid under the assumptions that the radiation 
is soft and that two-body collisions dominate the emission. It appears that conditions in 
SN1987a are such that radiation from the NN system is not soft enough to make the cal- 
culation we have presented here completely reliable. However, the framework described in 
this paper holds out the promise that corrections to our result can be systematically com- 
puted. In this way, it represents a significant advance on previous work, where uncontrolled 
approximations for the A^A^ dynamics were employed. Some of the higher-order corrections 
to the NLO result which was used to set the bounds (|5^ ) can still be expressed in terms of 
A^A^ system observables, e.g. derivatives of the A^A^ differential cross section with respect to 
angle and energy. Indeed, since the A^A^ differential cross section is essentially constant with 
energy and angle in the region probed by our calculation one might hope that the ©(x*^) 
contribution to £ is fairly small. Even were higher-order corrections to halve the emissivity 
the bound on the radius for the n = 2 case would still be a severe -R < 1 x 10^^ mm. This is 
two orders of magnitude beyond the bound attainable at any existing collider ||l3l. The only 
comparably stringent bound on the radius of toroidally-compactified, large, GODs comes 



from the cosmological considerations discussed by Hall and Smith in Ref. ||3^, who advocate 
i? < 5.1 X 10~^ mm. However, as Hall and Smith themselves comment, this bound is subject 
to their ignorance of the effect that GODs have on cosmology. 

In summary, for the case of two or three large, toroidally-compactified "gravity-only" 
dimensions one of the best bounds on modifications in the physics of gravity comes from 
delving into the extreme conditions present in SN1987a, and asking what effects such changes 
would have there. Our detailed examination of the two-nucleon dynamics which yields this 
bound implies that the radius of these extra dimensions must be below 7.1 x 10^^ mm for 
n = 2 and below 8.5 x 10^'' mm for n = 3. 
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